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Abstract 

Stochastic quantization is applied to derivation of the equations for the Wilson loops and 
generating functionals of the Wilson loops in the N = oo limit. These equations are treated 
both in the coordinate and momentum representations. In the first case the connection of the 
suggested approach with the problem of random closed contours and supersymmetric quantum 
mechanics is established, and the equation for the Quenched Master Field Wilson loop is derived. 
The regularized version of one of the obtained equations is presented and applied to derivation of 
the equation for the bilocal field correlator. The momentum loop dynamics is also investigated. 
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1. Introduction 



One of the possible ways of investigation of QCD is its reformulation in the space of all possible 
contours (the so-called loop space), when one considers gauge fields as chiral fields, defined on 
this space 1 ' 2 . Loop equations (for a review see 3 ) yield natural description of the confining phase 
of a gauge theory in terms of elementary excitations, which are nothing, but interacting closed 
strings 4 , and the main purpose of these equations is to choose in the large- N limit those of the 
free string theories, which describe gauge fields. This problem is not completely solved up to now, 
though several attempts to this direction are done (see for example 5 ). Loop equations are known 
to be nonlinear Schwinger equations in functional derivatives, formulated, generally speaking, for 
the generating functionals of the Wilson averages 6 , and that is why the solution of these equations 
is a difficult mathematical problem, which nowadays is disentangled only in two dimensions 7 
(notice also, that recently in 1+1 QCD the average of arbitrary number of Wilson loops on an 
arbitrary two-dimensional manifold was evaluated, and the relation to the string theory of maps 
was established 8 ). 

Recently alternative gauge-invariant equations were derived 9 using the stochastic quantiza- 
tion method 10 (for a review see 11 ), and investigated in various field theories 12 . These equations 
are written not for Wilson loops, but for field correlators and therefore are closely connected to 
the Method of Vacuum Correlators 13,14 (for a review see 15 ). The mathematical structure of the 
obtained equations is simpler than the structure of the loop equations, since the equations for 
correlators are just evolutionary integral-differential Volterra type-II ones. The derived system 
of equations for correlators is in some sense similar to the Bogolubov-Born-Green-Kirkwood-Ivon 
chain of equations for the Green functions in statistical mechanics 16 , connecting correlators with 
various number of fields, and therefore the complete information about the field theory under con- 
sideration is contained only in the full infinite set of equations. In order to obtain a closed system 
of equations one needs to exploit some kind of approximation. The irreducible field correlators (the 
so-called cumulants 9 > 12 > 13 > 14 >i5,i7 vanish, when any distance between two points, in which the fields 
in the cumulant are defined, increases and therefore it occurs natural to cut the obtained infinite 
chain of equations, using as a parameter of the approximation the number of fields in the cumulant 
(see discussion in 9 ). For example, in 9,12 the so-called bilocal approximation 13 ' 14 ' 15 , when all the 
cumulants higher than quadratic were neglected (which corresponded to the Gaussian stochastic 
ensemble of fields 17 ), was exploited in order to investigate the obtained equations perturbatively 
both in the standard, non-gauge- invariant, and gauge-invariant ways, perform stochastic regulari- 
zation, separate perturbative gluons' contributions, include external matter fields and apply the 
suggested approach to quantization of classical solutions. 

The aim of this paper is to apply the stochastic quantization method to derivation of equations 
for the Wilson loops in the N = oo limit alternative to the familiar loop equations. After that, 
using the cumulant expansion 9 ' 12 ' 13,14 ' 15,17 , one may get equations for correlators. 

There are two types of equations for the Wilson loops, we are going to present. The equation 
of the first type, which will be investigated in Sections 2 and 3, is the evolutionary equation 
of the heat transfer form with the functional Laplacian standing for the ordinary one. This 
equation is true in the limit N = oo for the arbitrary values of the Langevin time. It will be 
derived in Section 2 for a single non-averaged Wilson loop and then generalized to the case of 
the generating functional for the Wilson loops. After that we smear the functional Laplacian in 
the obtained equation, using the method, suggested in 18 , and reduce the problem to the integral 
Volterra type-II (by the Langevin time) equation. It will be shown that the action, over which 
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we shall perform the averaging in order to invert the smeared functional Laplacian (which is the 
action of the Euclidean harmonic oscillator at finite temperature 18 ), after performing the limiting 
procedure, which restores the reparametrization invariance, has a Gaussian form, and therefore 
for the problem under consideration the Langevin equation, where the role of the Langevin time 
plays the parameter of a contour, may be written. Thus we come to the problem of random 
closed paths 2 , which due to the fact that it may be described by the Langevin equation with 
the Langevin time being the proper time in the ensemble of contours is possible to be reduced 11 
to the supersymmetric quantum mechanical problem 19 , so that the dynamics of this system is 
governed by the two Fokker-Planck Hamiltonians 20 ' 21 , both of which survive (in contrast to the 
case considered in 20 ), since one should not take the limit of the Langevin time tending to infinity. 
Also the equation for the Wilson loop of the Quenched Master Field 22 is derived in Section 2, and 
it is discussed that due to the discontinuities of the Wilson loop in the momentum representation 23 
only the quenching prescription in the Langevin time direction yields nontrivial equation. 

In Section 3 we rewrite the equation for the non-averaged Wilson loop in the momentum 
representation, using the methods suggested in 23 ' 24 . We show that the role of the extra proper 
time introduced in 23 , in which the Wilson loop propagates and splits, is played by the Langevin 
time. 

In Section 4 we investigate the equation of the second type, which is also true only in the 
N = oo limit, but, in contrast to the equation of the first type, only in the asymptotical regime, 
when the Langevin time tends to infinity. This is the equation for the averaged Wilson loops, 
where the averaging is performed over the stochastic Gaussian noise fields, which is known to be 
equivalent at the Langevin time tending to infinity to the averaging with the gluodynamics action 
weight. However, it should be mentioned that since we work in the N = oo limit and neglect all 
the insertions of the stochastic Gaussian noise fields into the Wilson loop, this equation does not 
reproduce planar graphs. Up to now we do not know an equation for the averaged Wilson loops 
alternative to the familiar loop one, which reproduce correctly the planar graphs in the large-iV 
limit, where the functional Laplacian may be replaced by the operator of differentiation by the 
Langevin time. 

The equation obtained is then also generalized to the case of the generating functional 6 for 
the averaged Wilson loops, which provides Veneziano topological expansion 25 , and rewritten in 
the momentum representation. Then we regularize the equation in the coordinate representation, 
using the method suggested in 26 , after which it occurs to be defined on the space of smooth closed 
contours, to which one can apply the nonabelian Stokes theorem and the cumulant expansion. 
While integrated over the Langevin time t, this ultraviolet-regularized equation with the properly 
chosen t-dependent momentum cut-off yields the expression for the one-gluon-exchange diagram, 
where the perturbative gluon's propagator is written in the Feynman-Schwinger path integral 
representation 14 , and the Langevin time plays the role of the Schwinger proper time. Finally we 
use the obtained regularized equation in order to derive the equation for the bilocal correlator, 
whose Kronecker structure is supposed to be given. 

The main results of the paper are summarized in the Conclusion. 



2. Investigation of the heat transfer type equation in the coordinate 
representation 



We shall start with the Langevin equation 



10,11 
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K = (V A F AM )° - igrfr, (1) 
where V M = + [A^, ■} is the adjoint covariant derivative, F^ v = d^A v — d u A^ + [A^,A V ] is a 
strength tensor of the gluonic field, g = where A is the bare coupling, which remains finite 
while performing the expansion, and rf is a stochastic Gaussian noise field. 
Let us introduce the following loop functional of the Stokes type 

*n(C,t) = ^trT Pexpjdz^ (A„(z,t)+ig J dt'^t')] = 

c V o / 

= ^ tr A Ji^ IT ^1 + Atj ' s ^2 l {/S.Zj) l j i (2A lx (zj,ti)8{ti — t)+ igr]n(zj,ti))^ . 

Here in addition to the usual P-ordering we have introduced T-ordering because of the second 
term in the exponent, which contains t- integration. 

Differentiating ^ V (C, t) by the Langevin time one gets by virtue of (1) the following equation 

*„(C,t) = jjtrfdx^V x F^(x))T Pexpjdz^ (A^(z,t)+ig J dt'rj^t') 

C C V i 

= —trjdx^ x F x ^x))U{C,t)+igjdz^{C X0Z ,t) J dt' V ^z,t')<P{C ZX(V t)- 



c c 



-g 2 jdz^j dy v </>(C X0Z ,t) j dt' % {z,t')<p{C zy ,t) J dt" Vu {y,t")<f ) {C yXQ ,t) + ...], (2) 

CCO ' 

where (p(C xy ,t) = P exp j dz^A^z^t), and xo is an arbitrary but fixed point belonging to the 

Cxy 

contour C. 

It was proved in 6 that the expansion corresponded to the WKB approximation around the 
"classical" field in the effective scalar theory in the loop space. We see that this is really the case 
according to equation (2) since in the stochastic quantization method 10 ' 11 the degree of quantum 
correction of the grand ensemble of fields A^ and 77^ is determined through the maximal number 
of the noise fields entering the physical quantity under consideration, and on the right-hand side 
of equation (2) any new term of the expansion gives rise to the additional power of the noise 
field. 

In what follows we shall derive and investigate the equation for the pure "classical" field in 
the loop space ^(C, t) = jjtr<fi(C,t) corresponding to N = 00. This equation reads 

tt(C,t) = <[ dx,d*—^--y(C,t), t < +oc, (3) 

OCT\^{X) 



c 



where d x x {a) = J da' ' 



Notice that equation (3) is true for an arbitrary t < +00, and from the mathematical point 
of view it is a homogeneous functional heat transfer equation (which means that the loop space 



4 



Laplacian stands there for the ordinary one), while the loop equation in the N = oo limit is an 
inhomogeneous functional Laplace one. 

Using equation (3) one can easily write an equation for the generating functional of the Wilson 
loops 

Z(J(C)) = exp(^J(C)^(C) ] j, 

where the sum over the loops is defined as follows 6 

+00 

£/(<?)= If / Dx(a)f(x(a)). 

c +0 x(0)=x(T) 

This equation reads 

dtJJ{C) ~ I dX ^J^)JW) + ^ J[C ) t dz ^5a^(z)5J(C)6J(C>) (4) 

and does not depend on N explicitly. We shall see in Section 4 that this dependence appears 
in the equation for the generating functional of the averaged Wilson loops, which describes an 
"open string" with quarks at the ends. This equation will correspond to the Veneziano topological 
expansion ' . 

Let us now perform the smearing of the functional Laplacian, standing on the right-hand side 
of equation (3), by making use of the method suggested in 18 . It may be done via performing 
the polygon discretization procedure and then going to the continuum limit or directly using the 
continuum smearing without any reference to the polygon discretization. This smearing implies 
that one replaces the functional Laplacian 

1 cr+O 



A = f dxnffi- — — — da da'- — - — -- — — 

J ocf\^{x) J J_ o dx^(a') dx^(a) 



with the smeared one 



A (G) = [da v.p. f da'G(a - a') - - 6 + A, (5) 
J J Sx^a') 5x„(a) 

where G(a — a') is a smearing function, and the first term on the right-hand side of equation (5), 
which involves the principal-value integral, is an operator of the second order (does not satisfy the 
Leibnitz rule in contrast to the operator A) and reparametrization noninvariant. The averaging 
over the loops £(<r) is defined as follows 

< F[$ >f>= «^ , (6) 

e(o)=«i) 

where 
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, 1 1 

S = - J da J da' (a*)G- l (a - a')^')) , (7) 



and G 1 is an inverse operator. For the simplest case 

G(a - a') = e - , e < 1 
the action (7) reduces to the action of the Euclidean harmonic oscillator at finite temperature 
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S = \jda (et 2 {o) + -/{*)) , (8) 
o 

which may be recovered from the discretized action, when the number of vertices of the polygon, 
which approximates the loop, tends to infinity. The contribution of the first term on the right-hand 
side of equation (5) is of the order e for smooth contours since the region \a — a'\ ~ e is essential in 
the integral over da'. Hence when e tends to zero, this term vanishes, reparametrization invariance 
restores, and — > A. 

Thus equation (3) takes the form 

A iG) ^[x,t] = 4>[x,t], (9) 
where x = x(a) is the position vector of the contour C. Using the equation of motion 



< >f>= J da'G(a - a') < 



5F[£] (G) 



which follows from the oscillator action (8), one can solve equation (9) 18 with the condition 
\I/[0, t] = 1. Integrating over the Langevin time we arrive at the following equation 



J dt'^x, t'\ + l - J d 1 (<n x + vre, *] >f ) - < nv^t a >f 



+oo 

= t + - J d 1 (<^[x + ^,0] >f - < *[V7C, 0] >f } ) , (10) 
o 

which is the integral Volterra type-II one by the Langevin time, and the dependence on the initial 
conditions is put to the right-hand side. 

Finally let us go to the limit e — > in order to restore the reparametrization invariance of 
equation (10). Then performing the rescaling ( = we see that the averaging in both sides 

i 

of equation (10) becomes Gaussian with the action S = \ J daQ 2 (a). It means that we come to 

o 

the problem of motion of random closed paths, which occurs to be described by the Langevin 
equation, where the role of the Langevin time plays the parameter of the contour (i.e. the proper 
time in the loop space). The problem of random closed contours was investigated in 2 , where it 
was shown that the number of contours of the length T was given by the formula 

4N(T) = fHrf, 
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where D was the dimension of the space-time, c was some constant, and e was a cut-off parameter. 

Within our approach we reduced this problem to the Langevin equation, written in the proper 
time of the system under consideration. Therefore it is equivalent 11 to the supersymmetric quan- 
tum mechanical problem 19 . Namely if one writes down in components the equation of the Brownian 
motion of a classical particle in a heat bath, which describes the evolution of the position vector 
of a contour, 

dx 5S 

77" = + ^ 
oa ox 

then making the change of variables r] — > x in the partition function for the Langevin dynamics, 
we obtain the following effective Lagrangian 

where " ■ " = ' " = ^, V = S'. In order to have a closed supersymmetry algebra, one 

needs to introduce an auxiliary field D: 

so that the action, corresponding to (11), is invariant under the supersymmetry transformations 

Sx = sifj — ipe, 5D = eip + ipe, 5ip = (x + D)e, Sip — e(x — D). 

The dynamics of the system may be shown 20 ' 21 to be governed by the two Fokker-Planck Hamil- 
tonians 

h±- 1 62 +V 2 +V 

H --2Sx^ + 8 V ± 4^' 
where the Hamiltonian H~ corresponds to the propagation forward in the Langevin time a, and H + 
corresponds to the backward propagation. The existence of the two Fokker-Planck Hamiltonians 
is a consequence of the supersymmetry. However it should be emphasized that while in the 
usual causal interpretation of the Langevin equation, which requires propagation forward in the 
Langevin time, in the equilibrium limit only the forward dynamics survives since only the zero 
ground state of H~ contributes, and H + has strictly positive eigenvalues 20 , in our model it is not 
the case, because one should not take the limit a — > +oo. Therefore the dynamics of the system 
is governed by both the forward and backward Hamiltonians. 

To conclude this Section, we obtain the equation for the Quenched Master Field 22 Wilson loops. 
The important property of the Wilson loops in the momentum representation is the presence of 
finite discontinuities induced by the emission and absorption of gluons 23 

Ap^Si) = p^Si + 0) - p^Si - 0), 
and therefore only quenching in the Langevin time direction is valid 

Af(x, t) = e i(p5a - pM) 'A^(x), r]f{x, t) = e l(p5a ~ P5b)t f]^ b {x), 

where 
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< % b (x)f)?(y) >= 2^-S bc 5 a %J(x - y). 



Introducing a new loop functional 

Sl(C,t) = jjtr j dx^P lB^x,t)exp j> dy x A x (y,t)\ , 



c \ c 



where B® b (x,t) = i{p<~, a ~ Vbb)A ab (x,t), we obtain from (3), performing the smearing of the func- 
tional Laplacian as it has been done above, the following equation for the loop functionals \1/ and 

+00 

t] = 1 - - / d 1 (< si[ x + t\ >f> - < t] >f>) . 



3. Heat transfer type equation in the momentum representation 

In this Section we shall rewrite equation (3) in the momentum representation in order to 
investigate free momentum loop dynamics in the Langevin time. During this free motion the 
loop behaves as a collection of free particles, while the loops' collisions, at which momenta are 
rearranged between particles, which means that the loops split, will be considered in the next 
Section. 

The loop functional in the momentum representation is defined as follows 23,24 : 



= J DCexp ^ J p^dx^ *(C,t). 

Using the representation 



1 si+O ^ 2 

A = j ds\ j ds2 



for the functional Laplacian and following the procedures suggested in 23 ' 24 , one gets from equation 
(3) two equations, which describe the free propagation of the momentum loop: 

np,t) = -j:(^(s l )) 2 np,t), (12) 

i 

1 si 

*(P,t) = -2 J d Sl p'^ Sl ) J ds 2 p' tM (s 2 MP,t). 


Notice that in equation (12) the reparametrization invariance is preserved, i.e. when the parameter 
s changes to f(s), where f'(s) > 0, the positions of the discontinuities Sj shift, but their order and 
the values of Ap^Si) remain the same. 

Therefore one may conclude that the role of the extra proper time H, which was introduced 
in 23 , in which the Wilson loop propagates and bits, plays the Langevin time t. This evolution is 
however alternative to the evolution of the loop in the space of random contours, described by the 
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supersymmetric quantum mechanics (see the previous Section), where the role of the proper time 
played the parameter of the contour. 



4. An equation for the averaged Wilson loops 

This Section is devoted to investigation of the equation alternative to equation (3), which will 
be written for the averaged Wilson loops. It can be derived by averaging both sides of equation 
(3) over the stochastic noise fields 77" which in the physical limit of the Langevin time tending to 
infinity is equivalent 10 ' 11 to the averaging with the gluodynamics action weight. Therefore in the 
asymptotical regime t — > +00 one may use the equation of motion for the Wilson loop ty(C,t), 
which is nothing, but the loop equation in the limit N = 00. Thus we obtain the following 
equation: 

$(C,t) = xjdx^ v.p. J dy^8{x-y)^(C xy ,t)^(C yx ,t), t -> +00, (13) 
c c 

where $(C, t) =< ^(C,t — > +00) >^a, and the principal value integral on the right-hand side 
of equation (13) implies that we integrate only over those y^ 's, which are another points of the 
contour C as the point x. In other words the right-hand side of equation (13) does not vanish 
only when the contour C has self-intersections. 

As it was already mentioned in the Introduction, this equation does not reproduce the planar 
graphs since we have neglected insertions of the stochastic noise fields into the Wilson loop in the 
limit N = 00. It should be understood as an asymptotical at t tending to infinity equation for the 
"classical" Wilson average. 

Equation (13) may be generalized to the case of the generating functional Z(C,j) = ^ ^jj0r- , 

which is also a functional of the source j(C) = '^jp-- If N tends to infinity, but the ratio p = ^- is 
fixed, which corresponds to the Veneziano topological expansion 25 , generalizing the expansion, 
then the source j remains finite. The equation for the functional Z reads 

^Z(C) = kjdx, v.p. J dy„6(x - y) (n Z(C, y )Z(C.,,) + ^^jgy) + 

This is the Schwinger equation in functional derivatives, which describes the evolution of the wave 
functional of the open string with quarks at the ends (see discussion in 6 ). 

Let us now rewrite equation (13) in the momentum representation. It may be done using the 
procedure, suggested in 23 , and the answer has the form 

$(P,t) = -2A / ds 2 -—- [ d Sl -—-^pW,mP^,t), (14) 
where and are the parts of P from Si to s 2 and from s 2 to s x respectively: 

PM : p CQ( s ) = _ S21 ) s)j < s < s', p£\s) = - s 21 )s + S 2 i), s' < s < 1, (15) 
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P( 2 ) ; pW( 8 ) = Plx ( Sl + S2l5 ), < 8 < 1, (16) 

where s 2 i = s 2 - s u s' = 

Equation (14) describes the vertex of the splitting of the loop into two loops P — > P^ + p( 2 \ so 
that the momentum p^(s) of the initial loop is distributed between these two new loops according 
to the formulae (15) and (16). As was discussed in 23 , at every splitting the effective number of 
degrees of freedom per loop diminishes, so that if one approximates the initial loop by a polygon, 
the number of vertices of each of the polygons, approximating the loops after the splitting, will 
be smaller. At each splitting of the loop the function p M (s) smothers, and the derivative p'^s) 
decreases. Finally the loop, propagating in the Langevin time and splitting during this process, 
runs out of degrees of freedom and reduces to the point in the momentum space Pn(s) = p ii = const. 

Our next aim is to apply to equation (13) the nonabelian Stokes theorem and the cumulant 
expansion. After that using the bilocal approximation 13 ' 14 ' 15 , we shall derive the equation for 
the bilocal correlator, whose Kronecker structure will be supposed to be given. To this end we 
perform an ultraviolet regularization of equation (13), exploiting the method, suggested in 26 , 
where this procedure was done with the help of the heat-kernel regularized Langevin equation. 
The regularized version of equation (13) has the form 



y 

$(C, t) = A <j> <j> dy^ J Drexp 



c c 



( 1 A" 2 \ 

-- / dsf 2 (s) 
4 J 

V o / 



$(C xy r yx ,t)$(C yx r xy ,t), (17) 



where r(0) = x, r(A 2 ) = y, and the momentum cut-off A may be taken for example to be of the 
order of the inverse correlation length of the vacuum 15 T g ~ 0.2/m : |A| ~ However in what 

follows we shall use as a momentum cut-off 

It is known 2 that in the stochastic quantization of gauge theories it is not necessary to add 
a gauge-fixing term into the Langevin equation, since direct iterations in powers of the coupling 
constant without introducing ghost fields yield the same results as the Faddeev- Popov perturbation 
theory, because the Langevin time t takes the role of a gauge parameter. If one fixes t, calculates 
gauge-invariant quantities and then goes to the physical limit t — > +oo, all the lineary divergent 
in t terms will cancel each other in the same way as the terms, depending on the gauge parameter 
in the framework of the usual approach. 

Therefore our choise of the momentum cut-off looks natural from the point of view of the 
general principa of the stochastic quantization method. If we now integrate equation (17) by the 
Langevin time, analitically continuing both sides down from the asymptotics t — > +oo, then the 
obtained equation yields the expression for the second order one-gluon-exchange diagram in the 
perturbation theory in the nonperturbative gluodynamics vacuum 14 , where gluonic propagator is 
written in the Feynman-Schwinger path integral representation, and the role of the Schwinger 
proper time plays the Langevin time. This perturbative gluon propagates along a regulator path 
r xy , dividing the surface, swept by the nonperturbative gluodynamics string, into two pieces. 

Equation (17) is defined on the space of closed smooth contours, and therefore one can apply 
to it the nonabelian Stokes theorem and the cumulant expansion, after which the integral over 
the regulator paths, standing on the right-hand side of this equation, yields the amplitude of 
propagation of a free gluon. Let us vary both sides of equation (17), to which we applied the 
cumulant expansion in the bilocal approximation, by the element of the surface, lieing on the 
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contour C, and suppose for simplicity that the bilocal correlator has the only one Kronecker 
structure: 

< F flu (x,t)(f)(T xy ,t)Fxp(y,t)(/)(Ty X ,t) >= jj{8^\8 vp - 8 w 8 vX )D((x - y) 2 , t), 

where Y xy and Y yx are some contours. Then we get the following equation for the correlation 
function D: 

/Q f f 6 4* 

do Xp {x 2 )D{{ Xl -x 2 )\t) = — 1 -^jdx ll f dy»j—y 2 , 

s c c 

where the point x\ belongs to the contour C, the surface S is bounded by this contour, and the 

initial condition D(x 2 ,0) = is implied. 

5. Conclusion 

In this paper we applied the stochastic quantization method to derivation of the equations for 
Wilson loops in the N = oo limit. There are two types of equations: first is the heat transfer type 
equation (3), which is true for arbitrary values of the Langevin time. This is the evolutionary 
equation for the non-averaged Wilson loop. It was derived in Section 2 and then generalized to the 
case of the generating functional for the Wilson loops. After that, applying the smearing procedure 
to the functional Laplacian which stands on the right-hand side of this equation, we reduced the 
problem to the Volterra type-II integral equation (by the Langevin time) and established its 
connection with the problem of random closed paths. The latter occured to be described by the 
supersymmetric quantum mechanics of a particle in a heat bath, where the role of the proper time 
played the parameter of a contour. Therefore the dynamics of the system of random contours is 
governed by the forward and backward Fokker-Planck Hamiltonians. 

The quenching prescription may be applied to the obtained equation only in the Langevin time 
direction, since in the momentum representation the Wilson loop possesses discontinuities. The 
equation for the Wilson loop of the Quenched Master Field is presented at the end of Section 2. 

In Section 3 the heat transfer type equation was investigated in the momentum representation, 
and two alternative equations, describing the free propagation of the non-averaged Wilson loop in 
the Langevin time, were derived. 

In Section 4 an alternative equation for the averaged Wilson loops, which was also true in the 
limit N = oo but in the asymptotical regime of the large values of the Langevin time t, was derived 
and generalized to the case of the generating functional for the Wilson averages. This equation for 
the averaged Wilson loops, while rewritten in the momentum representation, describes splitting of 
the loop during its propagation in the Langevin time. At these splittings the loop loses its degrees 
of freedom and finally shrinks to a point in the momentum space. 

Then in order to apply to the equation in the coordinate representation the nonabelian Stokes 
theorem and the cumulant expansion and derive the equation for the bilocal cumulant, we regu- 
larized this equation, using the method suggested in 26 , after which, choosing for the momentum 
cut-off and integrating over t, one may recognize in the obtained expression the one-gluon- 
exchange diagram 14 , contributing to the averaged Wilson loop. Applying to the both sides of 
this equation the nonabelian Stokes theorem and the cumulant expansion and assuming that the 
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bilocal correlator had a given Kronecker structure, we obtained the equation for the corresponding 
coefficient function. 
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